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Abstract 

The observed luminosity distributions of X-ray sources indicate the presence of several populations of X- 
ray binaries in the nearby galaxies. Each population has its formation and evolutionary history, depending 
on the host environment. The features seen in the log N(>S) - log S curves for different types of galaxies 
and for different galactic components can be reproduced by a birth-death model, in which the lifespans of the 
binaries are inversely proportional to their X-ray brightness. Conversely, the dynamical history of a galaxy can 
be inferred from the luminosity distributions of its X-ray binary populations. 

1 Introduction 

X-ray binaries are powered by a compact star, which may be a neutron star or a black hole, accreting mate- 
rial from its companion. Systems having a massive OB companion star are called high-mass X-ray binaries 
(HMXBs), and systems with a low-mass companion star are known as low-mass X-ray binaries (LMXBs). In 
HMXBs the compact star accretes via capturing the stellar wind from its companion; mass transfer in LMXBs 
occurs when the companion star overflows its Roche lobe. 

X-ray binaries in the active states are luminous and easily detected within the Galaxy. Many X-ray binaries 
were also found in nearby galaxies by the Einstein and ROSAT X-iay satellites (Fig. 1, see also e.g. Fabbiano 
1995; Roberts & Warwick 2000), and recently, more are discovered by the Chandra X-ray observatory (see 
Weisskopf et al. 2000 for the description of the observatory). It is now evident that galaxies similar to our own 
(e.g. M31, Supper et al. 1997) normally host hundreds of active X-ray binaries. The current sample of X-ray 
binaries in external galaxies is sufficiently large that not only is reliable population analysis for an individual 
galaxy possible but we can also study their formation and evolution in different galactic environments. 

In this paper we present a simple birth-death model and calculate the populations of X-ray binaries in 
external galaxies. We also use the model to identify the relevant processes that give rise to the features in the 
luminosity distributions of X-ray binaries (such as the log N(>S) - log S curves) in different types of galaxies 
and in different galactic components. Here, only the basic formulation is presented and a few simple cases are 
shown as an illustration. Results of a more comprehensive study will be reported in Wu et al. (2001). 



2 X-ray binaries in a galaxy 
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ROSAT PSPC Observations of M81/M82 
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Figure 1: ROSAT PSPC image of the spiral galaxy M81 and the starburst galaxy The biightest point sources in M81 can be easily distinguished 

from the background. 



2.1 Accretion luminosity 

The luminosity of an accreting compact star is given by 
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where G is the gravitational constant, M and R are the mass and radius of the compact star, and M is the mass 
accretion rate. Typically, neutron-star binaries have luminosities Lx ~ 10'^^ erg s~^, but the luminosities of 
binaries with accreting white dwarfs would not greatly exceed lO'^'^ erg s~^. 

The expression for above is not strictly applicable to black-hole binaries, as black holes do not have 
a solid stellar surface that defines R. However, as accreting material attain speeds close to the light speed, c, 
when crossing the black-hole event horizon, we may parametrise the accreting luminosity as 
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where r/bh is the efficiency. If most of the X-rays are emitted from the inner accretion disk near the last stable 
orbit, then 7/bh ~ 0.1 for a Schwarzschild black hole and rybh ~ 0.4 for a maximally rotating Kerr black hole. 

The characteristic luminosity of neutron-star binaries can also be expressed as « rjas Mc^, because of 
the tight mass distribution of neutron stars. If we assume a canonical mass of 1.5 Mq for the neutron stars, then 
the value of the efficiency parameter ij^s is about 0.3. Hence, we can adopt ^ rj Mc^ (with 77 ~ 0.1) for 
both neutron-star and black-hole binaries, provided that the mass-transfer rates M are not significantly above 
the Eddington limit. 



2.2 Eddington limit 

The infall of the material onto a compact star is opposed by the accretion-generated radiative pressure, and 
the accretion process is self-regulating. The limiting luminosity Led, known as the Eddington luminosity, for 



spherical accretion is given by 

^--4.^— -i^j =1.3x10- ^—j ergs-, (3) 

where m,p is the proton mass, and ctt is the Thomson cross section. The Eddington luminosity depends linearly 
on the mass of the accreting star and is independent of other parameters. 

As the average luminosity of an accretor would not be greater than the Eddington luminosity, the lower Umit 
to the mass of the accretor can be constrained by the observed luminosity. Thus, we can use this as a working 
criterion to identify black-hole candidates. For instance, when we observe a source with X-ray luminosities 
greatly exceeding 2 x 10^* erg s~^, the Eddington luminosity of a l.S-M© compact star, we may classify it as 
a (candidate) black-hole system (if it is powered by accretion). 



2.3 Evolution and expected lifespan 

The duration of the active phase of an X-ray binary is limited by the lifespan of the companion star, and it is 
almost independent of the nature of the compact star. Massive OB stars have Ufespans of ~ 10^ — 10^ yr, 
implying that HMXBs powered by capturing the stellar wind from massive OB stars will not have a lifespan 
significantly longer than a few million years. The presently active HMXBs must therefore have been formed in 
very recent epochs. 

Low-mass stars have longer lifespans and evolutionary timescales than high-mass stars, and so LMXBs can 
remain X-ray active over a period longer than 10 yr. Mass transfer in LMXBs is usually driven by Roche-lobe 
overflow, which is caused either by the expansion of their companion stars when they evolve towards a red-giant 
phase, or by orbital shrinkage, when the binaries lose orbital angular momentum. 

For the first type of LMXBs, their mass-transfer rate is roughly given by (M) ~ M2/Tev, where M2 is the 
mass of the companion star, and t^v is the evolutionary timescale. The nuclear-evolution timescale of a star 
depends on its mass and evolutionary stage. The lifespan of F main-sequence stars is about 5 x 10^ yr; stars later 
than G type take more than 10^° yr to evolve away from the main-sequence stage. Clearly, nuclear evolution of 
main-sequence stars of masses ~ IMq is unable to sustain a persistent accretion luminosity ^ 10^^ erg s~^. 
Such accretion luminosities are possible for LMXBs with evolved companion stars. The duration of the red- 
giant phase is ^ 10^ yr for low mass stars (de Loore & Doom 1992), and hence. X-ray binaries with a red-giant 
companion are short-lived. However, galactic stars evolve into the red-giant stage at all epochs. X-ray binaries 
with red-giant companions are therefore formed continually throughout the life time of a galaxy. 

For the second type of LMXBs, orbital angular momentum can be extracted from the systems by emitting 
gravitational radiation (Paczynski 1967) or by rotational braking of the companion star via a magnetic steUar 
wind (Verbunt & Zwaan 1981). The timescale of orbital evolution driven by gravitational radiation is 

. ^ 5 c5 (Mi + M2)V3 
32(27r)8/3 G5/3 M1M2 

= 3^8xl0'.(l+,)n-(|^)""'(^)"%r ,4, 

(Landau & Lifshitz 1971, see also Wu 1997), where P is the orbital period. Mi and M2 are the masses of the 
compact star and the companion star respectively, and q (= M2/M1) is the mass ratio. As the timescale of 
orbital shrinkage is long (» 10^ yr), the corresponding mass-transfer rate is not large enough to account for 
the observed luminosity of the X-ray binaries. The orbital evolutionary timescale for magnetic braking is 

2 x IP'S fpy^ Me^T Ml pi°/3 
- (27r)W3 fc2 V (^^+^^)i/3M47^ec (5) 
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(Verbunt & Zwaan 1981; Wu 1997). The parameters k and / depend on the stellar model (see Skumanich 
1972), and the parameter 7 is determined by the mass-radius relation (R2 oc M^, where R2 is the radius of the 
companion). For low-mass stars, fc^ « 0.1, / ^ 1 and 7 w 1. The mass-transfer rate is higher than that in the 
case of gravitational radiation, and the active phase of these X-ray binaries usuaUy lasts ^ 10^ yr. 
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3 Populations of X-ray binaries 



In a simplistic point of view the populations of X-ray binaries in a galaxy observed at an epoch t are determined 
by the birth rates and the lifespans of the binaries. Assuming other factors are unimportant (for example, the 
two-body and three-body capture process in globular clusters, see Johnston & Verbunt 1996), we can construct 
a birth-death model and calculate the populations and luminosity distributions of X-ray binaries in a galaxy or 
in a particular galactic component. 

The basic formulation is as follows. Let the number of the X-ray sources at the luminosity range (L, L + dL) 
in a galaxy at time t be n{L, t)dL and the characteristic lifespan of these sources be r. The evolution of an 
X-ray binary population is governed by 

jn{L,t) = -kn{L,t) + f{L,t), (7) 

where k n{L, T) is the death rate, the f{L, t) is the birth rate, and (= r) is the characteristic lifespan of 
the sources. The number of sources with a luminosity brighter than L is simply 

/oo 
dL n{L, t) . (8) 

When the functional forms of k and /(L, t) are specified. Equations (7) and (8) can be solved easily. 
3.1 Impulsive birth 

X-ray binaries are bom (i.e. become X-ray active) when mass transfer from the companion to the compact 
star starts; and they die when the X-ray active phase ends. Suppose that there are two channels to create X-ray 
binaries in a galaxy: (i) systems that are formed continuously and the mass transfer of these systems is sustained 
by nuclear evolution of the companion star or by the orbital evolution of the binary, and (ii) systems that are 
born in a recent star-formation episode. The birth rate function f{L, t) thus consists of a steady (continual 
birth) component and an impulsive (starburst) component: 

= f,{L)[l + a6{t-t,)] , (9) 

where fo{L) is the birth rate of the steady component, a is the ratio of the strength of the impulsive component 
to the strength of the steady component, S{t — tg) is the Dirac 5-function, and ta is the time at which the 
starburst occurred. 

The effective duration of the active phase of an X-ray binary is determined by the mass of the companion 
star and the mass-transfer rate, and it is approximately given by M2/ (M). As M, the characteristic 

lifespan of an X-ray binary fc~^ can be parametrised as A; = (3L, where /3~^ = fj{M2)c^ and (M2) is the mean 
mass of the companion star, and 77 is a parameter which has the same order of magnitude as the parameters Tjbh 
and r]ns in §2. The fact that the characteristic lifespan of binaries with luminosities in the range {L, L + dL) is 
not expUcitly time-dependent allows us to integrate Equation (7) directly, yielding 



n{L, t) = no{L) e-^^* + ^ (^1 - e''^^*) + afo{L)e{t - ta)e 
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where 6{t — ta) is the Heaviside unit step function, and no{L) = n{L, 0) is the initial population. 

If the initial population and the birth rate are power-laws of L, i.e. rio = no*{L/L^)~"^ and fo = 
fo*{L/L^,)~°'^ (with power-law indices ai and 0:2 respectively), then 

N{>L,t) = no*L^{pL^t)°''-^T{l^ai,l3Lt) 



+ ^-j{l3L,tr- 
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+ afo,L4(3L4t - ta)r-^ r(l - a2, PL{t - ta))6{t - ta) , (11) 

where is a lower cut-off luminosity, and T{a, x) is the incomplete gamma function. In the limits of fo* 
(i.e. all the sources are primordial) the X-ray binary population has a distribution similar to the Schechter 
(1976) analytic luminosity function for galaxies in clusters. The same distribution is also obtained when 
{fo*/(3Lit) » rio* and a/3X* ^ 1 (i.e. most of the binaries were born in a recent starburst). The differ- 
ence between the Schechter (1976) luminosity function and the distribution that we obtain here is that in the 
latter case the luminosity break is time dependent — it is caused by the aging of the source population when it 
is not replenished. 
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Figure 2! Luminosity distributions of X-ray binaries in model galaxies for the impulsive-birth model. In the calculations we assume an age of 15 Gyr 
for the Universe. We also define two parameters Ac = (/o*/n-o*/3i») and Xa = {afo* /"-o*) to specify the relative strength of the continuous and the 



3.2 Cyclic generation 

In reality, star formation may not be a single event. For example, when two galaxies are orbiting around each 
other, periodic starbursts are induced by orbital interaction. In this case, the birth term /(L, t) has a periodic 
(starburst) component superimposed on a steady background: 

- ,UL)[l + b cos Luf] , (12) 

where 27r/a; is the characteristic timescale of the cycle. It follows that the population of X-ray binaries at time 
t is 
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and the number of systems brighter than L is 
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The integral I{a; u)t, (5Lt) above is defined as 
l{a;ujt,j3Lt) = I dx 
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(15) 
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and it has exact analytic form for integer a. The periodic component may introduce ripple-like features in the 
luminosity function, hence makes the log N(>S) - log S curves deviate from a simple or a broken power law. 



3.3 Other issues 

3.3.1 Neutron-star binaries 

The model above has not accounted for the fact that the persistent luminosity of neutron-star binaries cannot 
exceed the Eddington limit of a 1.5-Mq accretor. The luminosities of black-hole binaries are less restrictive, as 
there is no practical mass limit for black holes. A more proper treatment requires the populations of neutron- 
star and black-hole binaries be calculated separately. Although Equation (7) is applicable for both classes of 
objects, the corresponding expressions for N{> L) are not identical. For a population of X-ray binaries with a 
cutoff luminosity at L^d, N{> L) is given by 

N{> L,t) = J dLn{L,t). (16) 

The tight mass distribution of neutron stars together with the presence of a luminosity Umit will produce a spike 
at Led in the (differential) luminosity function n(L), and hence a cut-off in the observed log N(>S) - log S 
curve. 

The strength of the spike in n{L) is time-dependent. If all the neutron-star binaries in the galaxy were born 
at the same episode, because the systems with high mass-transfer rates cease to be active first and then those 
with lower mass-transfer rates, the spike will gradually decrease, on a timescale 

~ l.SxlO^^yr, (18) 

where Med and M^s is the Eddington mass-accretion rate and the mass of neutron stars respectively. For galax- 
ies with star-formation activity occurred more than 10^ yr ago, the spike has already degraded so substantially 
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that it may not be detected easily (Wu et al. 2001). The spike is, however, prominent for galaxies with violent 
starburst activity in the near past lO'^ yr), or if the neutron-star binaries are formed continuously. 

3.3.2 Transient sources 

The total duration of the X-ray active phases of transients is only some fraction of their hfespans. The parameter 
k is therefore not exactly equal to M/M2, but is proportional to it. As a first approximation, we may use an 
effective factor ^ to specify the fractional time at which the source is X-ray active, i.e. k = /3L/£_. If ^ does not 
depend on L explicitly, the solution to Equation (7) preserves its form. Otherwise, the solution and hence the 
luminosity functions are modified. 

It is worth noting that in addition to different fractional time ^ for their X-ray active phases, the transients 
and persistent sources have different discovery probabihties. This could complicate the calculations of the 
source populations and luminosity functions. Here, we have omitted the transient effects for simpUcity, but one 
must bear in mind that the effects might be important. 

4 X-ray sources in external galaxies 

Since 1999, Chandra has obtained high-resolution images of a number of nearby galaxies. These include 
the elhpticals M84 (NGC 4374) (Finoguenov & Jones 2001) NGC 4967 (Sarazin, Irwin & Bregman 2000), 
and Cen A (Kraft et al. 2000), the lenticulars NGC 1553 (Blanton, Sarazin & Irwin 2001) and NGC 1291 
(Irwin, Bregman & Sarazin 2001), the spirals M31 (NGC 224) (Garcia et al. 2000), M81 (NGC 3031) (Tennant 
et al. 2001), MlOl (NGC 5457) (Pence et al. 2001) and the Circinus galaxy (Sambruna et al. 2001), and 
the irregular starburst galaxy M82 (NGC 3034) (Matsumoto et al. 2001), and the mergers NGC 4038/4039 
(Fabbiano, Zezas (& Murray 2001). The Chandra observations show more than a hundred bright X-ray sources 
(with Lx > 10'^^ erg s~^) in the field each of these galaxies. While some of the bright sources are foreground 
stars and background AGN, the majority are probably X-ray binaries within the galaxies. The observations also 
indicate that the brightest sources in the early-type galaxies are generally brighter than their counterparts in the 
spiral galaxies, if the galactic nuclei are excluded. 

4.1 log N(>S) - log S curve 

The log N(>S) - log S curves of the Einstein and ROSAT sources in spiral galaxies are in general adequately 
fitted by a power law (see Fabbiano 1995). However, for the galaxies with a sufficiently large number of sources 
discovered (^ a few tens), the log N(>S) - log S curves are found to deviate significantly from a single power 
law, e.g. M31 (see Supper et al. 1997) and MlOl (Fig. 3, left panel; see also Wang, Immler & Pietsch 1999). 
A broken power law is often required in order to fit the log N(>S) - log S curves. The break is located at 
~ a few X 1037 erg s'^ (see e.g. M31, Supper et al. 1997; Shirey et al. 2001). 

The Chandra observations confirm the presence of the broken power-law in the log N(>S) - log S curves 
of X-ray sources in spiral galaxies (e.g. M81, Tennant et al. 2001). They also verify that the broken power 
law is present in the log N(>S) - log S curves for early-type galaxies, e.g. NGC 4697 (Sarazin, Irwin & Breg- 
man 2000), NGC 1553 (Fig. 3, right panel) and M84 (A. Finoguenov, private communication). The fact that 
a luminosity break is found in both the cases of early- and late-type galaxies suggests that the break may be 
universal. A possible explanation is that the luminosity break is caused by a population of neutron-star bina- 
ries which have super-Eddington mass-transfer rates. These neutron-binaries would have X-ray luminosities 
roughly about 2 x 10^^ erg s~^, the Eddington luminosity of the I.5-M0 accretors. If the location of the lu- 
minosity break is the same for all galaxies, it can be used as a distance indicator (Sarazin, Irwin & Bregman 
2000). 

However, when the X-ray sources in the disk and in the bulge of M81 are considered separately, very 
different log N(>S) - log S curves are obtained (Fig. 4, Tennant et al. 2001). The luminosity break is obvious 
in the log N(>S) - log S curve of the bulge sources, but the log N(>S) - log S curve of the disk sources is a 
single power law. Interestingly, the power-law slope of the disk sources is similar to that of the bulge sources 
below the luminosity break. 
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Ml 01 ROSAT HRI sources 





Figure 3: (Left) The log N(>S) - log S curve of the ROSATHRI sources in MlOl (NGC 5457). The data are taken from Table 3 in Wang, Immler & 
Pietsch (1999). The unit of S is counts per kilo-second. For an assumed distance of 7.7 Mpc to the galaxy, 1 unit corresponds roughly to a luminosity of 
2.8 X 10^* erg at 0.5-2.0 keV (1 cts ks'^ is approximately 4.0 x 10"^"' erg s=l). (Right) The log N(>S) - log S curve of the Chandra sources in 
NGC 1553. The data are taken from Table 1 in Blanton, Sarazin & Irwin (2000). The unit of S is counts per 10 kilo-second. The conversion between count 
rate and the unabsorbed X-ray luminosity is 4.91 x 10^^ erg cts~^ at the 0.3-10 keV (see Blanton, Sarazin & Irwin 2000 for details). 
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Figure 4: The log n(>s) - log s curves of the disk and bulge sources of M81 observed by Chandra. The sources are background subtracted with the 
log N(>S) - log S distribution of the background taken from Chandra calibration observation of CRSS J0030. 5+2618. Adopted from Tennant et al. (2001). 

4.2 Dynamical evolution of the host galaxy 

The absence of the break in the log N(>S) - log S curve of the disk sources in M81 is inconsistent with 
the proposition that the luminosity break is universal for all galaxies. The locations of the breaks at different 
luminosities for early- and late-type galaxies (cf. M81 and NGC 4697, see Tennant et al. 2001 and Sarazin, 
Irwin & Bregman 2000 respectively) suggest that the natures of the breaks are different for galaxies with 
different Hubble types. We propose that for some galaxies the luminosity break is determined by the formation 
and evolutionary history of the populations of X-ray binaries in the galaxy. The different morphologies for the 
log N(>S) - log S curves for the disk and bulge sources in the spiral galaxies reflect the different star-formation 
histories for the galactic components and the different evolutionary paths of their X-ray binary populations. 

We interpret the absence of the luminosity break in the log N(>S) - log S curve for disk sources in the 
spiral galaxy M81 as being a consequence of a continuous, smooth star-formation process. In normal spiral 
galaxies, star formation is triggered by compression of gas when the density wave sweeps across the galactic 
disk. The star-formation process is continuous and is relatively smooth (in comparison with the environment 
in starburst galaxies and mergers). As the X-ray binaries are formed continually in the spiral arms, the bright 
sources are dominated by young HMXBs with substantial proportion of black-hole binaries. The disk sources 
may also consist of two types of older LMXBs: one with their companion stars evolving into the red-giant 
phase, another one with their orbital separation sufficiently shrunk that the systems become semi-detached. 
Many of these LMXBs may contain accreting neutron stars. However, if the young HMXBs dominate the total 
population, there will not be a prominent luminosity break in the log N(>S) - log S curve. 

Formation of X-ray binaries in a galactic component is continuous if the galaxy is undisturbed. In this case, 
the population is dominated by LMXBs which are formed via processes the same as those of the LMXBs in the 
galactic disk. If the bulge is disturbed, e.g. by tidal interaction with another galaxy, a starburst is triggered by 
a sudden infall of gas, and a new population of stars as well as X-ray binaries are formed. As shown in §3.1, 
such an impulsive formation can produce a prominent luminosity break in the log N(>S) - log S curve, and the 
location of the break. When the population of the X-ray binaries ages, the break migrates to low luminosities, 
thus marking the look-back time of the catastrophic star-forming episode. 

Interestingly, the luminosity break in the log N(>S) - log S curve for the sources in the ellipticals (e.g. 
NGC 4697) seems to occur at a higher luminosity than that for the bulge sources in the spirals (e.g. M31 and 
M81) (A. F. Tennant and D. A. Swartz, private communication). If luminosity breaks are created solely by the 
aging of X-ray binary populations, the break seen in the log N(>S) - log S curves of the Chandra sources in the 
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early-type galaxies indicate that these galaxies had experienced some, possibly violent, star-formation activity 
in the very near past. 

However, this interpretation must be taken with caution. As pointed out in §3.3. 1, the tight mass distribution 
of neutron stars together with the presence of an Eddington Umit for each accretor could produce a spike in n{L) 
and hence a luminosity cut off in A^(> L). A substantial proportion of neutron-star binaries in the X-ray binary 
population can also cause a luminosity break in the log N(>S) - log S curve. The break caused by a population 
of neutron-star binaries with massive companions is expected to be located at ix ~ 2 x 10^^ erg s~^, the 
Eddington luminosity of a I.S-Mq accretor. Its strength will decrease on a timescale ^ 10^ yr, if there is no 
replenishment of the systems. As the luminosity breaks seen in the log N(>S) - log S curves of the sources in 
the bulge of the spirals scatter over a range of luminosities, they are unlikely to be all due to a population of 
accreting neutron stars. It is, however, possible that the breaks in the log N(>S) - log S curves of the early-type 
galaxies are caused by a population of neutron-star binaries with high mass-transfer rates. 

4.3 Super-Eddington sources 

Although super-Eddington sources seem to be more abundant in the early-type galaxies, given the different 
stellar contents in the galaxies, the excess of bright sources in the early-type galaxies is yet to be confirmed. 
The finding that the brightest sources in the spirals tend to be located in the disk (e.g. in M81, Tennant et al. 
2001) impUes that the disk sources are probably young systems. This also implies that the brightest sources 
in the spiral galaxies and the brightest sources in the early-type galaxies may not form in the same way. The 
population of bright sources in the arms of the spiral galaxies are dominated by HMXBs, in which the mass 
donors are young massive OB stars with strong winds; while the population of bright sources in the ellipticals 
and lenticulars are LMXBs, with a red-giant companion overflowing its Roche lobe. 

Despite the conventional scenario that X-ray sources in galaxies are mostly X-ray binaries, the super- 
Eddington sources may be an inhomogeneous class of exotic objects: hypemovae, supemovae, supersoft 
sources, intermediate-mass black holes, large clusters of massive OB stars, .... etc. These sources have not 
been included in our calculations presented above, and hence the results that we obtain may not truly reflect 
the real populations of X-ray sources in the galaxies. However, in order to take into account these objects, one 
must have some working models for their formation and evolution, which are unfortunately not available at 
present. On the bright side, if we accept the simple model in §3 by faith, we can explain the basic features in 
the log N(>S) - log S curves, and also use the observed populations of X-ray sources to probe the dynamical 
history of their host galaxy in the near past. 

5 Summary 

We have constructed a simple birth-death model and calculated the luminosity function of X-ray binaries in 
external galaxies. The model reproduces the features, such as the luminosity break, in the log N(>S) - log S 
curves of spiral galaxies. By assuming a continual formation process, the model also explains the absence of a 
luminosity break for the disk sources in M81. The location of the luminosity break depends on the look-back 
time of the previous starburst/star-formation activity, and so it can be used to probe the dynamic history of the 
host environment of the X-ray sources. 
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